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Abstract 

The current analysis of atmospheric neutrinos by the Super-Kamiokande Col- 
laboration is based only on fully-contained one-ring events and partially contained 
events. We show that the up-down ratio of fully-contained, inclusive, multi-ring 
events gives an independent test of the atmospheric neutrino anomaly, without the 
need for particle identification. Moreover, this class of events is rich in neutral 
current events and hence gives crucial information for discriminating between oscil- 
lations of into fe,T and Vs- 
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1 Introduction 



Recent results from the Super-Kamiokande collaboration have confirmed that the mea- 
sured neutrino fluxes, produced by cosmic ray showers in the Earth's atmosphere, do not 
agree with expectations from theoretical calculations. An anomaly is seen in events with 
both low visible energy, the sub-GeV data and high visible energy, the multi-GeV 
data 0. The collaboration has used three event topologies to present evidence for these 
anomalies: fully contained (FC) 1 ring events, where the Cerenkov ring has the charac- 
teristics of an electron or muon - e-like and /i-like events; and partially contained (PC) 
events, which Monte Carlo calculations show have a 98% probability to be produced by 
a Un charged current interaction. The e-like and /i-like events are studied both in the 
sub-GeV and multi-GeV samples, while the PC events are multi-GeV events. From the 
observed number of these events, the /^/e ratio, relative to Monte Carlo expectations, is 
measured to be 0.61 ± 0.03 ± 0.05 for sub-GeV data, and 0.66 ± 0.06 ± 0.08 for muhi-GeV 
data. The uncertainties are statistical and systematic, respectively. In the sub-GeV case 
the /i/e ratio is the ratio of 1-ring /x-like to 1-ring e-like events, while in the multi-GeV 
case the numerator also includes the PC events. 

The Super-Kamiokande collaboration has also measured the zenith angle distribution 
of these event classes, for both sub- and multi-GeV data. The e-like events were consistent 
with Monte Carlo expectations, while the distributions for /x-like and PC events were far 
from Monte Carlo expectations. Let Pi{4>) be the number of events of class i with visible 
products moving in the upward direction relative to the number of these events in the 
downward direction. The visible products must be moving in a cone about the vertical 
with half angle 0. These up-down ratios are important, as they are insensitive to much 
of the theoretical uncertainties associated with calculations of the neutrino fluxes from 
cosmic ray showers; however, they are only useful if there is a good angular correlation 
between the neutrino and the visible products of the neutrino interaction. For the multi- 
GeV data, the correlation between neutrino and charged lepton directions is good, with a 
mean of 15-20". For the sum of /i-like and PC events, Super-Kamiokande has measured 
p^(78°) = 0.52±0.07±0.01, (Monte Carlo expectation: 0.98±0.03±0.02), while from the 
e-like events they measure pe(78°) = 0.84 ± 0.13 ± 0.02 (Monte Carlo expectation: 1.01 ± 
0.06 ± 0.03), with uncertainties strongly dominated by statistics. While the measured 
value of pe is consistent with Monte Carlo calculations, the discrepancy in is especially 
significant. For the sub-GeV data, the correlation between neutrino and charged lepton 
directions is much poorer - with a mean of about 60° - and the data give p^(53°) = 
0.69 ± 0.08 and pe{53°) = 1.23 ± 0.12, where only statistical uncertainties are shown. 

This data, together with that of previous experiments, provides substantial evidence 
that neutrinos produced in the atmosphere are oscillating as they traverse the Earth. The 
low systematic uncertainty on p^, both from theoretical calculations of the flux and from 
the Super-Kamiokande detector, make it clear that the survival probability for a i/^ to 
traverse the Earth is substantially less than unity. We are therefore led to study other 
signals in the Super-Kamiokande data which could 
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• Confirm oscillations with low systematic uncertainties. 

Even though large uncertainties cancel from the /i/e ratio, several remain: the flux 
calculation, the charged current cross section, the neutral current cross section, the 
energy scale and the separation of 1-ring from multi-ring events each have system- 
atic uncertainties of 4-6% 0|. On the other hand, it is striking that the systematic 
uncertainties for the up-down ratios pe^^ are much smaller: the flux calculations, 
detector asymmetries and downward going muon flux each give systematic uncer- 
tainties of only 1-2% 0. The anomaly in the z/^/z/g flux ratio depends crucially on 
the particle identification, which has been quite convincingly tested at KEK 0] , but 
still one would like to have a test independent of the systematic issues of particle 
ID. 

• Probe the flavor of the oscillation mode. 

The measurement of clearly shows that are disappearing on traversing the 
Earth, and hence must have oscillated into a combination of u^, Vt and Vg (a singlet 
neutrino) .[] The present data strongly disfavors oscillations purely to i/g, but does 
not distinguish between pure Pr and Vg- A high statistics separation of the Pr and 
Vs modes requires event classes with substantial probabilities of being produced by 
neutral current interactions. 

In this letter we study the neutrino oscillation signal in the up-down ratio of inclusive 
multi-ring events, pMR, both in the sub- and multi-GeV data. These are all events in 
which neutrino interactions produce two or more charged particles or photon showers 
which lead to identified Cerenkov rings. The theoretical systematic uncertainties are 1- 
2% from flux calculations, and 1-4% from cross sections, and we expect the experimental 
systematic uncertainties also to be small. The large number of multi-ring events ensures 
small statistical errors; for example, in the multi-GeV data the number of multi-ring 
events is very similar to the sum of the number of 1-ring /x-like and PC events, so the 
statistical power of Pmt?** be the same as for p^. For the sub-GeV data the angular 
correlation will not be as strong as for the multi-GeV case, but the multi-ring events have 
a higher neutral current component, so that p^^i ^ay have better capabilities to separate 
v-r from Vs- This separation can also be probed via p^o, the up-down asymmetry of events 
with a single vr", as has been proposed While such events have a higher neutral 
current component, a measurement of p'm'r will not require any particle identification, 
and will have better statistics. It has been suggested that the and Vg modes could be 
distinguished by studying the fraction of 2-ring events which result from neutral current 
7r° production [Q. However, this method has large uncertainties due to the uncertainties 
in neutrino interaction cross sections. 

""^In this letter, we suppress our theoretical prejudice against singlet neutrinos, which apparently require 
either a new low energy scale of physics, or a non-minimal see-saw. 
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2 Up-down Ratios 



In this letter we present a simplified analysis for the up-down ratios pi{4>)- We assume 
that, for some choice of the cone half angle 0, the downward going neutrinos are essentially 
unoscillated, while the upward going neutrinos have completely oscillated.[| The measured 
zenith angle distributions make this a reasonable assumption. We also assume that the 
flux ratio of neutrinos produced in the atmosphere, r = z/^/z/g, is independent of energy, 
at least for the energies which dominate a given class of events.^ Then, for events of 
class i, which are induced by z/e charged current, charged current and neutral current 
interactions with relative probabilities /g, and f^c 

P^ = feiPee + rP^e) + U (P^^ + ^Pe^) + fNcil - Ps) (1) 

where Pij is the probability for oscillation z/j uj and Ps = {Pes + fP^s)/{^ + r). It is 
immediately apparent that, within these approximations, up-down asymmetries can only 
measure the combinations: Pee+rP^e, P^i^i + Pe^l/1^ and 1 — Ps; and the latter is only probed 
if /jvc is appreciable. The oscillation probabilities satisfy unitarity constraints, SjPjj = 1, 
but the three combinations which can be measured remain independent. The fractions 
/e, ffj, and fNc can be obtained from the Monte Carlo results of the Super-Kamiokande 
collaboration, and are shown in Table 0. 

In the cases of oscillations between two flavors, such as z/^ ^ z/g, z/^ — > z/^ or z/^ — ^ z/^, 
with CP conservation, the expressions for the ratio simplify drastically and involve only 
one parameter, P^^. For z/^ z/^ oscillations. Pee = 1, P/^e = Pe/x = -Ps = 0, and P^^ < 1, 
giving the up-down ratio 

P^^ = fe + UP,, + Inc. (2) 

For z/^ Us oscillation. Pee = 1, P,e = Pe, = 0, Pf,s = I - P,ti, and Ps = rP^s/i)- + r), 
and we find 

^ = /e + f,P,, + fNC^^. (3) 

]^ _|_ ^ 

Finally, the case of Ve is the most complicated one. Under the assumption of the 

CP conservation, P^g = Pe, = 1 — P,,-, and Pgg = P^^. We have P^ = in this case, and 
find 

pf = fe{r - (r - 1)P,,) + + ^^^^..) + fNC (4) 

In these cases, measurements of the up-down ratios for various event classes, i, should 
yield the same value for P^^. Values for P^^ extracted from p^"^'"^«'*« and p*^''-™"'** are 

^ Complete oscillation means that the relevant Am^L/E factors are sufhciently large that on averaging 
over L/E the probabilities become independent of these factors. If this is not the case, pi is still given 
by (0), with the probabilities now understood to be suitably averaged over energy. 

■^This is a better approximation for the sub-GeV data than for the multi-GeV case. When r has 
significant energy dependence, the up-down ratios are still given by (1), but with r and l/r suitably 
averaged. 
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fe 




Jnc 


#MC 


#data 


Fully Contained Sub-GeV (Analysis A) 


1 ring e-like 


0.88 


0.02 


0.10 


812.2 


983 


1 ring yU-like 


0.005 


0.96 


0.04 


1218.3 


900 


multi-ring 


0.24 


0.43 


0.33 


759.2 


696 


Fully Contained Multi-GeV 


1 ring e-like 


0.84 


0.07 


0.09 


182.7 


218 


1 ring /i-like 


0.005 


0.99 


0.005 


229.0 


176 


multi-ring 


0.30 


0.55 


0.15 


433.7 


398 


Partially Contained 


0.01 


0.98 


0.01 


287.7 


230 



Table 1: Data and Monte Carlo results reported for 25.5 kiloton-years, 414 days of running 
time, by the Super-Kamiokande collaboration |1|, 0] , using the flux calculations of Ref . . 
For each class of event, the first three columns give the Monte Carlo results for the 
fractions of the events which are induced by z/g charged current, charged current and 
neutral current interactions, the fourth column the expected number of events, and the 
fifth column the measured number of events. 

shown in Table 2, using r = 2.15 for sub-GeV data, and r = 3 for multi-GeV data 
For oscillations to z/^ or z/^, P^^ is essentially independent of r. For oscillations to z/g, as 
r is increased from 2.5 to 3.5, the value of P^^ varies from 0.20 to 0.32 (from p™"'*'), and 
from 1.13 to 1.08 (from p™"'**). There is good consistency for oscillations of z/^ to Ur or 
Vsi but oscillations to z/g are excluded by the multi-GeV data at greater than 6o", for any 
reasonable value for r. 

For oscillations of z/^ to z/,- or z/^, the sub-GeV data appears to give a somewhat larger 
value for P^^ than does the multi-GeV data. However, the measured value of p'*"''(53°) is 
larger than the true value, pt-, due to the smearing effect from the poor angular correlation 
of the sub-GeV data. This can be parameterized by a phenomenological parameter D: 
p*"*(53°) = pt + D{\ — Pt)- The central values of P^^ extracted from sub- and multi-GeV 
data coincide ii D 1/3. This is consistent with the quoted angular resolution for the 
sub-GeV data. 

To measure Pmri it is necessary to define a direction for the multi-ring events, and this 
can be done in many ways. In contrast to single-ring and PC events, it cannot be defined 
using one particular ring; the direction must be defined in an inclusive way. The vertex 
should first be determined using the standard method, then vectors from the vertex to 
individual hits in the photomultiplier tubes (PMTs) can be drawn. One way to define 
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l^fl l^T 


l^s 


l^fi l^e 


p-^(53°) 


0.67 ±0.08* 


0.68 ±0.08* 


0.38 ±0.15* 




Consistent at 2a* 


Consistent at 2cr* 


0.77 ±0.12* 




0.52 ±0.07 


0.52 ±0.07 


0.24±0.11 




Consistent at Icr 


Consistent at la 


1.10 ±0.08 



Table 2: Values of P^^ extracted from data on up-down ratios, assuming perfect angular 
correlation, includes both 1-ring /i-like and PC events. * Only statistical uncertainties 
are included in the sub-GeV case. 

the direction of the event is the sum of all the vectors weighted by the corrected number 
of photo-electrons in the PMTs An equivalent way is to use the analog of the thrust 
variable in the QCD jet studies at e^e~ colliders, 

T = max^— — -— — , (5) 

Ei \Pi\ 

where the maximum is obtained by varying a vector with a unit length n. The n which 
maximizes the thrust is the direction of a jet. At water Cerenkov detectors, the momentum 
vector. Pi, can be taken as the vector to each PMT weighted by the corrected number 
of photo-electrons in the PMT. This value is approximate, as it ignores the dependence 
of photo-electron production on particle type 0. The advantage of using the thrust 
variable is that one can study the thrust distribution to see if Monte Carlo calculations 
give a reasonable agreement with the data. Whatever definition of the direction of the 
event is employed, a detector simulation is necessary to study the correlation between the 
direction of the event and the original neutrino direction. 

3 Multi-GeV Data 

We propose an analysis of the zenith angle dependence in inclusive, multi-GeV, FC multi- 
ring events. According to the Monte Carlo study by the Super-Kamiokande Collaboration 
summarized in Table |l], the FC multi-ring events are dominated by CC events. Since 
the multi-ring event sample is richer in multi-pion production and hence to higher energy 
neutrinos than the single ring event sample, the angular resolution is expected to be 
good. We actually do not know the breakdown of the FC multi-GeV multi-ring event 
sample into quasi-elastic (QE), single pion, and multi-pion production; the corresponding 
breakdown in the sub-GeV sample for Kamiokande experiment by Kajita shows that 
the QE contribution is small (less than 5.9% with CC, CC and NC combined), 
and the rest is roughly equally divided between single pion and multi-pion production. 
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In the multi-GeV sample, the multi-pion production is expected to be more important, 
and can roughly be approximated by the Deep Inelastic Scattering (DIS) processes, even 
though the energy of the neutrino is still relatively low. For the DIS processes, the 
correlation between the neutrino direction and the momentum of the hadronic system 
and the charged lepton can be studied, and is of course perfect if one can observe all of 
the momenta. Even though one misses both particles below the Cerenkov threshold and 
neutrons, the correlation should be better than for QE events, where it is 15-20% (RMS) 
0. Therefore, we assume in this paper that the PC multi-GeV multi-ring events have 
reasonable angular resolution, at least as good as the QE events. We advocate evaluating 
p1Tr^{'P) "with a large half cone angle of = 78'', so that only about 20% of the events 
are not used. 

Once a direction of the event is defined, as outlined in the previous section, we can 
study the zenith-angle dependence of the PC multi-GeV multi-ring events. Since more 
than half of the events are expected to be z/^ CC events, the same anomaly seen in /x-like 
one-ring events must also appear in this zenith angle dependence. In a 25.5 kt-years data 
sample, 398 events were observed, which gives roughly 11% statistical uncertainty in the 
up-down ratio. 

The up-down ratio is also useful to study the nature of the oscillation, to distinguish 
T-'ii oscillation from oscillation. We denote the survival probability of z/^ 

as P^^ as in the previous section. In the case of — > oscillation, the rate for upward 
going CC events is suppressed by the factor P^^, while v^. CC and NC event rates are 
unchanged. The resulting up-down ratio can be calculated using Eq. (0) and Table |l|, 

Por Vg oscillation, Eq. gives 

1 + r 

Here, the v^^ to v^. flux ratio is roughly r ?a 3, suitably averaged over a range of energies 
and zenith angles, which can be studied in detail with their Monte Carlo analysis. Pinally, 
the case of — > v^, oscillation is more complicated. The up-down ratio is 

^(ew*. = 0.30(r - (r - 1)P^^) + 0.55 {^^ + ^^^P^/.) + 0.15 (8) 

from Eq. (H). 

One important question is the theoretical uncertainty in the composition of the PC 
multi-GeV multi-ring events. The thesis by Shunsuke Kasuga varied the CC and 
NC cross sections for sub-GeV data sample very conservatively by ±30% and ±50%, 
respectively. We take the same variation and determine the range allowed for the up- 
down ratio clS cl function of P^^. 
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Figure 1: Up-down ratios in the FC multi-GeV multi-ring events for three oscillation 
scenarios. The CC and NC cross sections are varied by ±30% and ±50%, respectively, 
to give bands for each cases. The data point with large error bars correctly reflects 
the observed up-down ratio in multi-GeV /i-like events and its error in the horizontal 
direction, while the vertical position is yet to be measured. The large vertical error bar 
is the expected statistical uncertainty using 414 days of data. The smaller error bars 
correspond to four times more data - about 5 years. 



Figure |^ shows the up-down ratios predicted for the three oscillation scenarios as a 
function of the survival probability P^^. An important point is that the up-down ratio 
is quite insensitive to the theoretical uncertainty in the cross sections. Another point is 
that three scenarios are relatively well separated, even though z/^ and Ug are somewhat 
close because of the small NC fraction. The current up-down ratio in the multi-GeV yU-like 
events determine P^^ with an accuracy shown by the large error bar. The up-down ratio 
in the FC multi-GeV multi-ring events is not reported; we simply took a number with an 
anticipated statistical error. It is clear that the current data sample can distinguish the 

case and the other two. The measurement of the up-down ratio in this data sample 
will provide an independent consistency check of the observed anomaly, without relying 
on any particle identification. 

4 Sub-GeV Data 

The sub-GeV multi-ring data is potentially even more interesting than the multi-GeV data 
because it is rich in neutral current (NC) events, with a fraction f^c of about a third. 
Therefore, this data sample has considerable power to discriminate between z/^ — > 
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oscillations and z/^ — > Vr oscillations. 

The angular resolution is expected to be poorer than for the multi-GeV data. However, 
judging from a Monte Carlo calculation of the zenith-angle distribution of the one- ring /x- 
like events with neutrino oscillations 0, the zenith angular dependence is not completely 
washed out; rather it is diluted only by D ~ 0.4. The Super-Kamiokande Collaboration 
quotes a mean angular correlation in one-ring events to be 54% for muons and 62% for 
electrons [|l|. Therefore, it appears that the zenith-angle dependence or, in particular, the 
up-down ratio can be studied even with the sub-GeV data sample. To be conservative, 
we study Pm5j(0) with = 53°, ie. we use only the up-most bin (cos 6 < —0.6) and the 
down-most bin (cos 6 > 0.6) of the data. In this case, good angular resolution is not 
needed: A(cosO) ^ 2/5, or A9 ^ 53° is sufficient, and is a reasonable expectation. 

The most interesting component of the sub-GeV, multi-ring events is the NC, which 
is made of single 7r° production and multi-pion production with almost an equal amount 
according to Kajita's table 0. The angular dependence in the single 7r° production was 
studied by Diwan and Goldhaber ||^, and is quite good. They proposed to study this mode 
exclusively by enhancing the angular resolution with a kinematical cut whose efficiency is 
lower than 30%. We, on the other hand, are interested in studying the multi-ring sample 
inclusively to collect maximum statistics. Even without cuts, the contamination from 
the wrong bins is not large. Furthermore, the multi-GeV /i-like data suggests that the 
zenith-angle distribution is more-or-less flat for the first two bins and the last two bins. 
This also helps us to reduce the problem of poor angular resolution. Finally, multi-pion 
production in the NC data is expected to show a good angular resolution, much better 
than the one-ring events. The NC multi-pion production in the sub-GeV category is 
coming from higher energy neutrinos than the single- ring events because (1) the loss of 
the lepton energy makes the NC events appear less energetic than they actually are, and 
(2) the cross section rises as a function of energy unlike the QE and single-pion processes. 
We studied the angular correlation between the neutrino momentum and the momentum 
of the hadronic system numerically, and found indeed that the correlation is excellent 
(Fig. In reality, one cannot measure the total momentum of the hadronic system, 
because neutrons are not seen, and protons have a high energy threshold for Cerenkov 
radiation. The true resolution needs to be studied with a full detector simulation. Still, 
it appears reasonable to assume that the angular resolution is sufficient for our analysis. 

Following the same analysis as for the multi-GeV data set, we show the expected 
behavior of the up-down ratio p'mr ^ function of the u^j, survival probability P^^ for 
three oscillation scenarios. By varying the CC and NC cross sections by ±30% and ±50%, 
respectively, the three scenarios give bands shown in the Figure 3. Just like for the multi- 
GeV sample, the ratios are quite insensitive to a large variation of the cross sections. The 
important point is the Ur and scenarios are well separated. This figure is drawn with a 
dilution factor D = 1/3 in the up-down ratio, conservatively allowing for a poor angular 
correlation. In this plot, we assumed that all of the events had the same dilution, which is 
probably a too pessimistic assumption; the multi-ring NC events must have better angular 
correlations than the QE events as discussed above. 
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Figure 2: The angular correlation between the momentum of the neutrino and the 
hadronic system W in DIS NC events. The neutrino energy spectrum is taken from 
the calculation of Honda et al . 



The published data set has 696 sub-GeV multi-ring events. Because we conservatively 
use only the up-most and down-most bins, the expected current statistical error is about 
12%. It is clear that oscillations to Ue and to can be clearly separated. In a few years, 
the and scenarios will be separated as well. 

5 Conclusions 

We have proposed a new test for atmospheric neutrino oscillations using currently available 
data: the up-down ratio of FC, multi-ring, inclusive events, pmr- This method has good 
statistics and does not require particle identification. We anticipate that, with 414 days 
of data, two signals for neutrino oscillation could be obtained, one at the 3a level. This 
would provide an important, independent test of atmospheric neutrino oscillations. These 
measurements, one with sub-GeV data and the other with multi-GeV data, could also 
confirm that z/^ oscillate to Ur or Ug, and not to z/g. We believe that Pmr, like pe,/x are 
currently statistics limited, so that the significance of these measurements will improve 
with time. Combining p^/ij ^^"^ Pma** results could allow a 2a separation of u^- and 
modes with 5 years of data. With 5 years data it is likely that Am^ can be reliably 
extracted from the zenith angle dependence of the 1-ring e-like, 1-ring /i-like and PC 
events. It may then be possible to increase the significance of the sub-GeV, multi-ring 
analysis by including data from all zenith angles. We believe that our sub-GeV analysis 
has been conservative: we have included a large dilution factor for angular correlations, 
we have taken very large uncertainties in the neutrino cross sections, and we have not 
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Figure 3: Up-down ratios in the FC sub-GeV multi-ring events for three oscillation sce- 
narios. The CC and NC cross sections are varied by ±30% and ±50%, respectively, to 
give bands for each cases. The data point with large error bars correctly reflects the ob- 
served up-down ratio in multi-GeV /^-like events and its error in the horizontal direction, 
while the vertical position is yet to be measured. The large vertical error bar is the ex- 
pected statistical uncertainty using 414 days of data. The smaller error bars correspond 
to four times more data - about 5 years. Only the upmost (cos O < —0.6) and downmost 
(cos0 > 0.6) bins out of 5 are used, and a dilution of the ratio is included by taking 
£>= 1/3. 
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attempted to optimize the cone angle 0, so that it will be possible to improve the power 
of this technique. 

Our analysis has largely been based on the assumption that the oscillations involve 
only two flavors of neutrinos. Ultimately, it will be important to perform a global analysis 
of the various up-down ratios using the predictions of (1) for general neutrino oscillations. 
The flt is overconstrained, as there are seven event classes: e-like 1-ring, /i-like 1-ring and 
multi-ring for the sub-GeV data; FC e-like 1-ring, /i-like 1-ring and multi-ring, as well as 
PC events, for the multi-GeV data. The averaged neutrino flux ratio, r, is close to 2 for 
the sub-GeV data, but for the multi-GeV data it is considerably higher, and depends on 
the cone angle 0. Hence, with several years of data, it will be possible to flt the seven 
up-down ratios to obtain values for the flve oscillation probabilities P^e, Pe/j,, P^e^ P^l^l and 

Ps. 
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